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A cubic kilometer scale experiment has been proposed to detect cosmic neutrinos of energy from 
tens of GeV up to the highest energies observed for cosmic rays, ~ 10 20 eV, or possibly even beyond. 
Detection efficiencies depend crucially on the neutrino-nucleon cross section at these energies at 
which radiative corrections beyond the lowest order approximation could become non-negligible. 
The differential cross sections can be modified by more than 50% in some regions of phase space. 
Estimates of corrections to the quantities most relevant for neutrino detection at these energies give, 
however, less dramatic effects: The average inelasticity in the outgoing lepton is increased from 
~ 0.19 to ~ 0.24. The inclusive cross section is reduced by roughly half a percent. The dominant 
uncertainty of the standard model ultra-high energy neutrino-nucleon cross section therefore still 
comes from uncertainties of the parton distributions in the nucleon at very low momentum fractions. 

PACS numbers: 13.15.+g, 13.60.Hb, 13.85.Tp, 95.55.Vj 
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I. INTRODUCTION 



Several proposals have recently been put forward for 
the search for cosmic neutrinos above tens of GeV up 
to the high energy end of the cosmic ray spectrum, and 
possibly beyond. The most well developed technique is 
to detect Cherenkov light from the muon produced in a 
charged-current reaction of an ultra-high energy (UHE) 
muon neutrino with a nucleon in either ice or water. Sev- 
eral prototype detectors based on this technique have 
been constructed, namely DUMAND at a depth of nearly 
5 km in the ocean near Hawaii (now out of commission) , 
the neutrino telescope at Lake Baikal about 1 km deep, 
NESTOR at about 3.5 km depth in the Mediterranean 
near Greece, and AMANDA up to 2 km deep in South 
Pole ice [0. Other techniques have been proposed such 
as detection of horizontal air showers or detection 

of acoustic (D or radio waves M associated with the neu- 
trino induced cascades. 

For a given neutrino flux, the efficiency of all these 
methods depends predominantly on the neutrino-nucleon 
interaction cross section. To lowest order in the elec- 
troweak (EW) coupling, this cross section has been dis- 
cussed in detail in the literature, see, e.g., Refs. ||,|7j for 
the most recent work. Due to uncertainties in the extrap- 
olation of quark distribution functions in the nucleon to 
very small fractional momentum transfers, x < 10 -4 , and 
large (negative) squared four-momentum transfers, Q > 
10 5 GeV 2 , best estimates for energies around 10 20 eV (in 
the laboratory frame) vary by factors of a few. 

On the other hand it is well known that higher or- 
der processes can become important or even dominant 



for electromagnetic (EM) interactions at UHE. For ex- 
ample, energy exchange between two leptons l\ and I2 is 
dominated by EM bremsstrahlung, h + h h + h + J, 
for center-of-mass energies s exceeding the square of the 
electron mass, rather than by ordinary Mott scattering, 
^1+^2 —* h + h- The relevant cross section rises with the 
logarithm of s. 

EW radiative corrections to deep inelastic neutrino- 
nucleon scattering have been calculated before in the 
literature: Ref. g] contains a discussion of the leading 
log approximation for which only corrections involving 
photons are relevant. These corrections contain a factor 
ln(s/m 2 ) where mi is the charged lepton mass and also 
depend on the behavior of the parton distribution func- 
tions. Single and double differential cross sections have 
been evaluated in Ref. || for neutrino energies E v up to 
a few hundred GeV, using rough estimates of the parton 
distributions available back then. 

To compare with HERA measurements Q at s ~ 
10 5 GeV 2 , corresponding to E v ~ 50TeV, more recently 
complete analytical expressions have been presented in 
Refs. 0,0. At these energies, corrections have been 
shown to be up to 50% in certain areas of phase space. 
Consequently, at energies approaching 10 20 eV, radiative 
corrections could be larger still and may play an impor- 
tant role. We therefore found it worthwhile to extend 
estimates of radiative corrections to UHE, using updated 
parton distribution functions and with a special empha- 
size on the quantities relevant for UHE neutrino detec- 
tion. 

The rest of this paper is organized as follows: In Sec. II 
we introduce our estimates of radiative corrections at 
UHE. In Sec. Ill we present numerical results. We fi- 
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nally summarize our findings and resulting consequences 
in Sec. IV. 



II. ESTIMATES OF RADIATIVE CORRECTIONS 
AT ULTRA-HIGH ENERGIES 

In principle, the full EW radiative corrections complete 
to order g 2 where g is the EW coupling constant can be 
computed from the expressions given in Ref. |1Q] and 
references therein. These, however, involve hundreds of 
terms and are somewhat hard to reproduce. In contrast, 
the leading logarithmic approximation is comparatively 
simple and is expected to be very good at the UHEs we 
are interested in Furthermore, to be concrete we will 
restrict ourselves to charged-current reactions which are 
more relevant for most of the detection methods relying 
on Cherenkov radiation from muons. We use the usual 
kinematic variables, Q 2 = 2ME u xy, with M the nucleon 
mass, y = (E v — Ei)/E u , and Ei the energy of the outgo- 
ing charged lepton in the laboratory frame. The contri- 
bution from the lepton leg to the double differential cross 
section in leading log approximation can then be written 
as 
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where a ~ 1/137 is the EM fine structure constant, 
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Zmm = l- y + xy, 



and d 2 ag/dxdy is the lowest order cross section. 
Eq. ([!]) is expected to be a good approximation for 
ln(2ME n u/m 2 ) » ln(l — y) Isj. Assuming isoscalar nu- 
cleons (i.e., averaging over protons and neutrons), the 
average of the lowest order cross section over neutrinos 
and antineutrinos is given by 
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Here, Gp is Fermi's constant, to^ is the mass of the W 
boson, and q v (Q 2 ,x) and q s (Q 2 ,x) are the valence and 
sea quark distributions of the nucleon. 

The first term under the integral in Eq. (Q) comes from 
bremsstrahlung of the lepton. Its z — ► 1 infrared diver- 
gence is cancelled by the second term under the integral 



which describes lepton propagator corrections due to vir- 
tual photons such that the total result is finite. It was 
argued in Ref. Q that the lepton contribution to the 
QED correction is dominant and that contributions from 
the quarks do not have to be included because they are 
accounted for in the quark distribution functions. 

For comparison we also calculated the hard photon 
bremsstrahlung contribution (Thb(E u ) to the inclusive 
cross section. In the following we denote the modu- 
lus of the 3-momenta of the outgoing photon and the 
quark by k and q' , respectively, and we also use q' Q — 
{{q') 2 + M 2 ) 1 / 2 . Furthermore, let p be the modulus of 
the 3-momentum of the incoming neutrino and fikp, Hkq', 
and Hp q i the cosine of the angle between the 3-momenta 
indicated. After integrating out the 5— functions for 4- 
momentum conservation we have 



<Jhb(E„) 



1 1 Mmi 

4(2^) 4 E 2 V 



J dkdq dfikpdfJ-kq' dx 



k 



xM~ 2 



1 2 qv(Q 2 ,x) + q s (Q 2 ,x) 



1/2 



(4) 



where fi pq > is given by 
q' { M 



k ( q' 



k , s M M(k~M)+mf/2 

+-1-/%,- — + — ri — > 5 

q q qxE y 

and the integration in Eq. (^) is performed over all ar- 
eas where — 1 < [i vq i < 1 and where the angle depen- 
dent factor in the integrand is real. The squared matrix 
element M 2 in Eq. (0) is averaged and summed over 
polarizations of incoming and outgoing particles, respec- 
tively, as well as over the electric charges of the quarks 
in the relevant distributions. In computing this matrix 
element we used the software REDUCE. The Feynman 
diagrams for bremsstrahlung photons attached to the in- 
coming and outgoing quark, the outgoing lepton and to 
the intermediate W boson as well as all interference terms 
were included. We note that as opposed to Eq. ([I]) this 
causes some double counting of radiative quark correc- 
tions, when using the same quark distribution functions 
as for Eq. ([!]). This is, however, expected to cause only 
a small error jll] . 

As expected, Eq. ([|) diverges logarithmically for k — » 
0. One can, however, still get a rough estimate of the 
radiative corrections without having to compute the vir- 
tual corrections. To this end one observes that a natural 
cutoff scale for k is given by the quark and lepton masses, 
and we will therefore restrict the range of integration to 
k > min(M, to/). This will not provide a precise value 
for the EM radiative corrections but at least a rough es- 
timate. 
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FIG. 1. The lowest order inclusive charged-current neu- 
trino-nucleon cross section as given by Eqs. (|^) and (^j) as 
a function of neutrino energy in the nucleon rest frame. 



For the quark distribution functions we use the CTEQ3 
distributions in the deep-inelastic scattering factorization 
scheme parametrization at Q 2 — (1.6 GeV) 2 Assum- 
ing a scaling of the sea quark distribution with InQ 2 , we 
adopt 
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This analytical form yields good agreement of the lowest 
order cross sections with Ref. |?]] (see Fig. Q). Although 
these distributions are uncertain within up to factors of 
a few for x <C 1CP 4 and Q 2 ;§> m 2 ^, our results for the 
fractional radiative corrections are relatively insensitive 
to these uncertainties. 

The cross sections are also practically independent of 
the lepton mass at these energies and therefore apply to 
all neutrino flavors although the resulting shower devel- 
opment is different for electrons and tau leptons. 



III. NUMERICAL RESULTS 

We now present results of numerical evaluations of 
Eqs. (|l|) and (||). We first note that for neutrino detection 
at UHE the probably most relevant quantity is the inclu- 
sive charged-current cross section cr(E u ) because, unlike 
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FIG. 2. The logarithmic distribution of y = 1 — Ei/E v 
from the lowest order calculation, Eq. (^) (dashed lines) and 
including radiative corrections calculated in leading logarith- 
mic approximation, Eq. (jl|) (solid lines) for E v = 10 20 eV 
(thick lines) and for E v = 10 5 eV (thin lines). 



in an accelerator experiment such as HERA, the pro- 
posed methods are not very sensitive to the phase space 
distribution of produced hadrons and charged leptons but 
rather act as calorimctric detection techniques. The only 
exception may be the quantity y which is the fraction 
of the neutrino energy initially going into the hadronic 
channel of the cascade ensuing, and equals 1 minus the 
fraction transferred to the outgoing charged lepton whose 
Cherenkov radiation is observed in the optical techniques. 
We therefore restrict ourselves to the quantities a{E v ) 
and y(da/dy){E v ,y)/(j{E v ). 

Fig. U presents the distribution of y, comparing the 
lowest order result from Eq. (||) and the result includ- 
ing radiative corrections in leading log approximation 
from Eq. (|l|) for a neutrino energy E v = 10 20 eV and 
E v = 10 5 eV. The resulting values for the average inelas- 
ticity (y) are 0.190 and 0.307 for the lowest order cross 
section, and 0.240 and 0.334 for the corrected one, re- 
spectively. We note in passing that corresponding values 
for the neutral-current neutrino-nucleon scattering cross 
sections are similar. The up to 3-dimensional integrals 
involved in these calculations have been evaluated using 
standard numerical techniques. 

Fig. |^ compares the fractional corrections 5{E V ) = 
a(E u )/ao(E u ) — 1 to the inclusive lowest order cross sec- 
tion resulting from Eqs. (|l|) and (Q) as functions of E v . 
To perform the 3- and 5-dimensional integrations, respec- 
tively, we used the adaptive Monte Carlo routine VE- 
GAS developed by Peter Lepage jl3| in the version from 
Ref. p4| . In order to ensure a reasonably smooth in- 
tegrand and sufficiently accurate results, we resorted to 
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FIG. 3. Absolute values of fractional radiative corrections 
8(E U ) = a(Eu) /(To {El,) — 1 to the inclusive lowest order cross 
section as functions of E v . The solid line shows the negative 
of the total correction in leading log approximation [Eq. (jl|)] 
and the dashed line is the hard bremsstrahlung contribution 
[Eq. (Q)] . Wiggles in the curves are due to the finite accuracy 
of the Monte Carlo integration which is at the few percent 
level. 



integration variables that are power laws in x and, typ- 
ically, logarithms in the other variables. Note that the 
leading log approximation for the total correction is neg- 
ative and roughly an order of magnitude smaller than 
the hard bremsstrahlung contribution which tends to be 
cancelled by virtual corrections. 



IV. DISCUSSION AND CONCLUSIONS 



As can be seen from Fig. 0, at energies around 
10 20 eV the radiative corrections to the single differ- 
ential charged-current neutrino-nucleon cross section 
(da/dy)(E l/ ,y) are negative for y < 0.02 and positive 
otherwise. For y < 10~ 3 they grow larger than 50%, 
whereas for y > 0.2 they are of the order of 30%. The 
average inelasticity (y) is increased from ~ 0.19 to ~ 0.24 
whose potential influence on the development of the neu- 
trino induced cascade is probably the strongest effect of 
radiative corrections. The corrections to the total cross 
section are negative and roughly constant at about half 
a percent (see Fig. |||). Apart from physics beyond the 
standard model (see, e.g., Refs 1 15 I3] , uncertainties of 
the UHE neutrino-nucleon cross section to date therefore 
are by an ample margin dominated by the uncertainties 
in the parton distributions in the nucleon at very low 
momentum fractions. 
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